COMGEN-BEM: Boundary element model generation for composite materials micromechanical analysis

Composite Model Generation-Boundary Element Method (COMGEN-BEM) is a program developed in PATRAN command language (PCL) which generates boundary element models of continuous fiber composites at the micromechanical (constituent) scale. Based on the entry of a few simple parameters such as fiber volume fraction and fiber diameter, the model geometry and boundary element model are generated. In addition, various mesh densities, material properties, fiber orientation angles, loads, and boundary conditions can be specified. The generated model can then be translated to a format consistent with a boundary element analysis code such as BEST-CMS

Composite micromechanical modeling using the boundary element method

The use of the boundary element method for analyzing composite micromechanical behavior is demonstrated. Stress-strain, heat conduction, and thermal expansion analyses are conducted using the boundary element computer code BEST-CMS, and the results obtained are compared to experimental observations, analytical calculations, and finite element analyses. For each of the analysis types, the boundary element results agree reasonably well with the results from the other methodologies, with explainable discrepancies. Overall, the boundary element method shows promise in providing an alternative method to analyze composite micromechanical behavior

An Improved Plastically Dilatant Unified Viscoplastic Constitutive Formulation for Multiscale Analysis of Polymer Matrix Composites Under High Strain Rate Loading

Polymer matrix composites are commonly used to fabricate energy-absorbing structures expected to experience impact loading. As such, a detailed understanding of the dynamic response of the constituent materials is necessary. Since the rate, temperature, and pressure dependence of carbon fiber reinforced polymer matrix composites are primarily manifestations of the rate, temperature, and pressure dependence of the polymer matrix, it is crucial that the constitutive behavior of the matrix be accurately characterized. In this work, an existing unified viscoplastic constitutive formulation is extended to ensure thermodynamic consistency and to more accurately account for the tension-compression asymmetry observed in the response of polymeric materials. A new plastic potential function is proposed, and elementary loading conditions are utilized to determine relations between model constants to ensure nonnegative plastic dissipation, a necessary thermodynamic requirement. Expressions for plastic Poissons ratios are derived and are bounded by enforcing nonnegative plastic dissipation. The model is calibrated against available experimental data from tests conducted over a range of strain rates, temperatures, and loading cases on a representative thermoset epoxy; good correlation between simulations and experimental data is obtained. Temperature rises due to the conversion of plastic work to heat are computed via the adiabatic heat energy equation. The viscoplastic polymer model is then used as a constitutive model in the generalized method of cells micromechanics theory to investigate the effects of matrix adiabatic heating on the high strain rate response of a unidirectional composite. The thermodynamic consistency of the model ensures plastic dissipation can only cause an increase in temperature. Simulation results indicate that significant thermal softening due to the conversion of plastic work to heat is observed in the composite for matrix dominated deformation modes

Utilization of the Generalized Method of Cells to Analyze the Deformation Response of Laminated Ceramic Matrix Composites

In order to practically utilize ceramic matrix composites in aircraft engine components, robust analysis tools are required that can simulate the material response in a computationally efficient manner. The MAC/GMC software developed at NASA Glenn Research Center, based on the Generalized Method of Cells micromechanics method, has the potential to meet this need. Utilizing MAC/GMC, the effective stiffness properties, proportional limit stress and ultimate strength can be predicted based on the properties and response of the individual constituents. In this paper, the effective stiffness and strength properties for a representative laminated ceramic matrix composite with a large diameter fiber are predicted for a variety of fiber orientation angles and laminate orientations. As part of the analytical study, methods to determine the in-situ stiffness and strength properties of the constituents required to appropriately simulate the effective composite response are developed. The stiffness properties of the representative composite have been adequately predicted for all of the fiber orientations and laminate configurations examined in this study. The proportional limit stresses and strains and ultimate stresses and strains were predicted with varying levels of accuracy, depending on the laminate orientation. However, for the cases where the predictions did not have the desired level of accuracy, the specific issues related to the micromechanics theory were identified which could lead to difficulties that were encountered that could be addressed in future work

Modeling of Damage Initiation and Progression in a SiC/SiC Woven Ceramic Matrix Composite

The goal of an ongoing project at NASA Glenn is to investigate the effects of the complex microstructure of a woven ceramic matrix composite and its variability on the effective properties and the durability of the material. Detailed analysis of these complex microstructures may provide clues for the material scientists who `design the material? or to structural analysts and designers who `design with the material? regarding damage initiation and damage propagation. A model material system, specifically a five-harness satin weave architecture CVI SiC/SiC composite composed of Sylramic-iBN fibers and a SiC matrix, has been analyzed. Specimens of the material were serially sectioned and polished to capture the detailed images of fiber tows, matrix and porosity. Open source analysis tools were used to isolate various constituents and finite elements models were then generated from simplified models of those images. Detailed finite element analyses were performed that examine how the variability in the local microstructure affected the macroscopic behavior as well as the local damage initiation and progression. Results indicate that the locations where damage initiated and propagated is linked to specific microstructural features

Associative Flow Rule Used to Include Hydrostatic Stress Effects in Analysis of Strain-Rate-Dependent Deformation of Polymer Matrix Composites

designing reliable composite engine cases that are lighter than the metal cases in current use. The types of polymer matrix composites that are likely to be used in such an application have a deformation response that is nonlinear and that varies with strain rate. The nonlinearity and the strain-rate dependence of the composite response are due primarily to the matrix constituent. Therefore, in developing material models to be used in the design of impact-resistant composite engine cases, the deformation of the polymer matrix must be correctly analyzed. However, unlike in metals, the nonlinear response of polymers depends on the hydrostatic stresses, which must be accounted for within an analytical model. By applying micromechanics techniques along with given fiber properties, one can also determine the effects of the hydrostatic stresses in the polymer on the overall composite deformation response. First efforts to account for the hydrostatic stress effects in the composite deformation applied purely empirical methods that relied on composite-level data. In later efforts, to allow polymer properties to be characterized solely on the basis of polymer data, researchers at the NASA Glenn Research Center developed equations to model the polymers that were based on a non-associative flow rule, and efforts to use these equations to simulate the deformation of representative polymer materials were reasonably successful. However, these equations were found to have difficulty in correctly analyzing the multiaxial stress states found in the polymer matrix constituent of a composite material. To correct these difficulties, and to allow for the accurate simulation of the nonlinear strain-rate-dependent deformation analysis of polymer matrix composites, in the efforts reported here Glenn researchers reformulated the polymer constitutive equations from basic principles using the concept of an associative flow rule. These revised equations were characterized and validated in an experimental program carried out through a university grant with the Ohio State University, wherein tensile and shear deformation data were obtained for a representative polymer for strain rates ranging from quasi-static to high rates of several hundred per second. Tensile deformation data also were obtained over a variety of strain rates and fiber orientation angles for a representative polymer matrix composite composed using the polymer

Finite Element Model for Failure Study of Two-Dimensional Triaxially Braided Composite

A new three-dimensional finite element model of two-dimensional triaxially braided composites is presented in this paper. This meso-scale modeling technique is used to examine and predict the deformation and damage observed in tests of straight sided specimens. A unit cell based approach is used to take into account the braiding architecture as well as the mechanical properties of the fiber tows, the matrix and the fiber tow-matrix interface. A 0 deg / plus or minus 60 deg. braiding configuration has been investigated by conducting static finite element analyses. Failure initiation and progressive degradation has been simulated in the fiber tows by use of the Hashin failure criteria and a damage evolution law. The fiber tow-matrix interface was modeled by using a cohesive zone approach to capture any fiber-matrix debonding. By comparing the analytical results to those obtained experimentally, the applicability of the developed model was assessed and the failure process was investigated

Micromechanical modeling of laminated composites with interfaces and woven composites using the boundary element method

The boundary element method is utilized to analyze the effects of fiber/matrix interfaces on the micromechanical behavior of laminated composites as well as the elastic behavior of woven composites. Effective composite properties are computed for laminated SiC/RBSN and SiC/Ti-15-3 composites, as well as a woven SiC/SiC composite. The properties calculated using the computerized tool BEST-CMS match the experimental results well

High Strain Rate Deformation Modeling of a Polymer Matrix Composite

Recently applications have exposed polymer matrix composite materials to very high strain rate loading conditions, requiring an ability to understand and predict the material behavior under these extreme conditions. In this second paper of a two part report, a three-dimensional composite micromechanical model is described which allows for the analysis of the rate dependent, nonlinear deformation response of a polymer matrix composite. Strain rate dependent inelastic constitutive equations utilized to model the deformation response of a polymer are implemented within the micromechanics method. The deformation response of two representative laminated carbon fiber reinforced composite materials with varying fiber orientation has been predicted using the described technique. The predicted results compare favorably to both experimental values and the response predicted by the Generalized Method of Cells, a well-established micromechanics analysis method

High Strain Rate Deformation Modeling of a Polymer Matrix Composite

Recently applications have exposed polymer matrix composite materials to very high strain rate loading conditions, requiring an ability to understand and predict the material behavior under these extreme conditions. In this first paper of a two part report, background information is presented, along with the constitutive equations which will be used to model the rate dependent nonlinear deformation response of the polymer matrix. Strain rate dependent inelastic constitutive models which were originally developed to model the viscoplastic deformation of metals have been adapted to model the nonlinear viscoelastic deformation of polymers. The modified equations were correlated by analyzing the tensile/ compressive response of both 977-2 toughened epoxy matrix and PEEK thermoplastic matrix over a variety of strain rates. For the cases examined, the modified constitutive equations appear to do an adequate job of modeling the polymer deformation response. A second follow-up paper will describe the implementation of the polymer deformation model into a composite micromechanical model, to allow for the modeling of the nonlinear, rate dependent deformation response of polymer matrix composites